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Abstract 

Study of Lya absorption lines in the spectra of 
double line of sight (DLOS) quasars holds the promise 
V) of diagnosing the nature of the structures that give 
rise to Lya absorption. Based on simulations of 
bJ[)DLOS with various separations through a single pop- 
^ ulation of absorbers, four tests have been designed to 
diagnose absorber mass distributions (smoothly vary- 
ing with radius or irregular), geometries (spherical or 
y—( disk/slab-like), and kinematics (isotropic or system- 
atic velocities). Applying the tests to existing data at 
z ~ 2 we find that: (1) The observed neutral hydrogen 
column density (Nh) distributions of Lya lines coin- 
cident to both LOS are consistent with a smooth mass 
distribution. (2) Observed large anticoincident Nh 



< 



are not consistent with a single population of smooth 
spherical absorbers, which should exhibit sharp cut- 
offs at small Nh in the Nh distribution. (3) There 
^j^is marginal evidence that the observed RMS value of 
I velocity differences between coincident lines increases 
O with DLOS separation as is expected for disk/slab- 
like absorbers that have systematic velocity fields. 
(4) The observed velocity dispersion along a single 
LOS is small compared to the RMS difference be- 
tween widely separated LOS, which is not consistent 
with models of cloudlets moving isotropically within a 
spherical structure. Overall, only a smooth disk/slab- 
like model with systematic velocities remains consis- 
tent with the inferred properties of a single population 
of Lya absorbers. 

Subject headings: quasars: absorption lines — quasars: 
general 
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Introduction and Motivation 

Observations of coincident and anticoincident ab- 
sorption lines 4 in the spectra of quasar pairs whose 
projected separations are small (hereafter double line 
of sight, DLOS) promise to place constraints on the 
nature of Lya absorbers (Sargent, Young, & Schnei- 
der 1982, Foltz et al. 1984, Smette et al. 1992, Bech- 
told et al. 1994, Dinshaw et al. 1994 (hereafter F84, 
S92, B94, and D94)). We focus on observations of the 
neutral hydrogen column density Nh and the velocity 
difference distributions of z ~ 2 Lya forest lines along 
DLOS quasars. Our goal in this Letter is to demon- 
strate that various geometric and kinematic models 
give rise to discernibly different DLOS observations, 
which provide diagnostics of the mass distributions 
and kinematics of the structures. We do not at- 
tempt detailed predictions based upon rigorous mod- 
els, but focus on idealistic, phenomenologically moti- 
vated models selected to illustrate expected behavior 
for more realistic scenarios. Monte-Carlo simulations 
of DLOS with various separations are presented for a 
spherical structure with an isotropic velocity field, a 
constant rotation "disk" , and a flattened "slab" with 
constant velocity flowing radially toward or away from 
a central point (planar inflow/outflow). The resulting 
model Nh distributions and velocity differences of the 
absorption lines are compared to current DLOS ob- 
servations. We argue that the most consistent single 
population model is a flattened disk/slab-like struc- 
ture with a Nh distribution that smoothly falls off 
with radius and with some level of systematic motion 
in the plane. 

The sizes of Lya absorbers, inferred from the frac- 
tion of lines with coincident redshifts in both spec- 
tra, are too large to be explained by most traditional 
models. Assuming a spherical population of identical 
structures, B94 and D94 derived a most probable ra- 
dius of i?~ 80 h' 1 kpc (h = #o/100 kms" 1 Mpc" 1 , 
q = 0.5) for Lya absorbers at z ~ 2 from their stud- 
ies of Q1343+266 A,B. At redshifts 0.5 < z < 0.9, 
coincident lines in the spectra of the more highly sep- 
arated (D ~ 340 h' 1 kpc) Q0107-025 A,B provide an 
even larger strict lower limit of D/2 on the absorber 
radius (Dinshaw et al. 1995). These low redshift data 
should be considered separately since it is a matter 



i Roughly, coincident lines are those which appear in both spec- 
tra of a close pair, whereas anticoincident lines are those which 
are observed in one of the spectra but have no detectable coun- 
terpart in the other spectrum. 



Charlton, Churchill, & Linder 



Probing Ly« Absorbers 2 



of current debate whether high and low redshift Lya 
absorbers comprise a single population. Since two of 
the four coincident lines in Q0107-025 A,B fall in the 
Lyman limit regime this letter will focus on the z ~ 2 
population. 

The distributions of equivalent widths and veloc- 
ity differences of coincident absorption lines at z ~ 2 
show evidence of becoming less correlated with in- 
creasing DLOS separation. UM673 A,B, a lensed sys- 
tem with LOS separations ranging from D ~ — 
2 h~ x kpc, yields k, 50 absorption lines in each spec- 
trum. All but two 5 of the absorption lines are coinci- 
dent along the two LOS, with equivalent widths and 
LOS velocities that are identical within errors (S92). 
Q2345+007 A,B has LOS separation D ~ 11 h' 1 kpc, 
assuming a lensed system with zi ens = 1.49 (Fischer 
et al. 1994). If we account for likely metal lines (Stei- 
del & Sargent 1991), eleven lines are coincident and 
one is anticoincident (F84). The equivalent widths 
of the coincident lines are less correlated than those 
of UM673 A,B. The LOS RMS velocity difference is 
25 km s _1 , but within the experimental errors this is 
consistent with km s" 1 (F84). Q1343+266 A,B, a 
physical pair with separation D ~ 40 kpc, yields 
?s 10 absorption lines in each LOS. Although the B94 
and D94 line lists differ in detail (likely due to dif- 
ferent spectral resolutions), significant variations in 
equivalent widths between the two LOS are observed 
(cf. Fig 3 of D94). D94 find an RMS velocity differ- 
ence of 65 ± 24 km s _1 among six coincident absorp- 
tion lines. 

In principle, coincident Lya absorption lines could 
arise either because a continuous structure covers 
both LOS, or because a "cluster" of smaller absorbers 
intercepts both LOS. Consider four idealized mod- 
els of Lya absorbers: 1) a continuous mass dis- 
tribution in a spherical structure, 2) discrete ab- 
sorbers ("cloudlets") distributed in a spherical struc- 
ture, 3) a continuous mass distribution in a disk/slab- 
like structure, and 4) discrete cloudlets distributed in 
a disk/slab-like structure. Based upon the large frac- 
tion of coincidences seen in Q1343+266 A,B (B94, 
D94), a cloudlet scenario requires near unit cover- 
ing factor for column densities greater than detection 
thresholds. For the disk geometry, a covering factor 
near unity implies a filling factor near unity, so that 
little difference exists between the latter two models. 

Continuous mass distributions likely result in a 



5 The two possible anticoincidences may be a Mgll doublet. 



monotonically decreasing Nn(r), where r is the radius 
of the structure. At the opposite extreme, discrete 
cloudlets within a structure could give rise to an irreg- 
ular Nh distribution, independent of r. The kinemat- 
ics within absorbing structures may range from highly 
systematic (e.g. rotation or planar inflow/outflow) to 
quite "random" (isotropic). To bracket a reason- 
able range we study two models with smoothly vary- 
ing Nn(r): (1) an isotropic velocity distribution in 
a spherical absorbing structure of radius R, which 
we interpret as discrete cloudlets moving isotropically 
with velocities ~ (GM/R) 1 ! 2 , where M is the total 
mass of the absorbing structure, and (2) disk/slab- 
like structures of radius R with either a constant 
rotation curve or a constant planar inflow/outflow 
with respect to a central point. We compare these 
models to a model with an irregular Nh distribu- 
tion. The smooth "isotropic sphere" serves to illus- 
trate mini-halo class models (Rees 1986, Miralda- 
Escude & Rees 1993), the "rotating disk" reflects 
the expected Lya absorption properties of Vanishing 
Cheshire Cat (VCC) galaxies (Salpeter 1993, Salpeter 
1995, Salpeter & Hoffman 1995), and the planar in- 
flow/outflow "slabs" loosely illustrate Zeldovich pan- 
cakes (Cen et al. 1994, Shandarin et al. 1995, Hern- 
quist, Weinberg, & Katz 1995). 

Monte Carlo Simulations and Model Results 

We assume all Lya forest absorption lines at z ~ 2 
are produced by a single population of absorbing 
structures. Thus, the observed column density dis- 
tribution f(N H )dN H oc N H P dN H must be recov- 
ered from the Nu(r) distribution across the absorb- 
ing structure, following Nu(r) oc r 2 ^ 1- ' 3 ). We use 
f3 = 1.5 (Petitjean et al. 1993, Rauch et al. 1993), 
which results in Njj(r) oc r -4 . Using Monte-Carlo 
techniques, we simulate absorption properties along 
two LOS (A and B). We use the parameter d = D/R 
to represent the ratio of the DLOS separation to the 
radius of the absorbing structure, and examine the 
values d = 0.01, 0.1, 0.25, 0.5, 1.0, and 1.5. For 
each LOS, we compute a normalized column density 
N = Nh(j)/Nh(R) and a normalized LOS velocity 
V = V\ os /V*, where Nh(R) is the neutral column 
density at the absorber "edge" , V\ os is a velocity vec- 
tor projected along the LOS, and V* is the velocity 
vector magnitude. For comparison with observations, 
Nh(R) represents the limiting column density of a 
sample. 
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Simulations for the disk geometry account for ran- 
dom orientations on the plane of the sky, which re- 
sult in elliptical cross-sections. In this case, Nn(r) 
is the column density at the radial distance on the 

1/2 

disk/slab, given by r = / |cos 2 a + (sin a/cos i) 2 j- , 
where / is the LOS distance from the ellipse center, 
i is the inclination subtended between the perpendic- 
ular axis of the disk and the plane of the sky, and 
a is the angle subtended between a vector pointing 
to the LOS and the semi-major axis of the ellipse. 
For rotational motion, the LOS velocity is given by 
V*(l/r) cos a sin i, where V* is the disk circular ve- 
locity. For planar inflow/outflow, V* represents the 
constant velocity magnitude and the LOS velocity is 
±14 (l/r) sin a tan i. 

In Fig. 1, we present the distribution log(A4i) ver- 
sus log(TVs), which we call /(NajNb)- An irreg- 
ular Nh distribution (independent of r), in which 
f(N A ,N B ) = N^^Ng 1 - 5 , is shown in panel (a). 
The smooth Nn(r) distribution is illustrated in pan- 
els (b-f), where panel (d) shows a spherical absorb- 
ing structure for d = 0.25 and panels (b), (c), (e), 
and (f) show randomly oriented disks for d = 0.01, 
0.1, 0.5, and 1.0. There is no expected correlation be- 
tween Na and Nb for the irregular Nh distribution, 
which yields an /(Na,Nb) independent of d. The 
smooth spherical structure yields an increased prob- 
ability density along the distribution "boundaries" . 
The higher density arises from the large range of ran- 
dom angles that locate LOS B at similar radii, yield- 
ing similar Nh- Random inclinations of the smooth 
disks result in (Na, Nb) pairs that lie outside the re- 
gions of increased probability density. For small d, 
the coincident LOS are at virtually the same radii 
within a structure, so that the Nh are highly corre- 
lated. As d increases, the correlation weakens until at 
d > 1 the second LOS no longer passes through the 
structure at Nb > Nh(R) if the first LOS has large 
Nh (small r). 

The Nh distributions of anticoincidences f(N ac ), 
shown in Fig. 2, exhibit remarkable differences be- 
tween the smooth spheres and smooth disks. For 
spherical structures with d < 1, the f(N ac ) sharply 
cut off at N ac = N H (R) (1 - d)~ 4 . For disks, the 
f(N ac ) follow a power law for large N ac , but exhibit 
increased probability densities at low columns that 
are more pronounced for smaller d. The latter arise 
due to increased probability for an anticoincidence in 
the case of intermediate to high inclination, when the 
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Fig. 1. — The distributions of column densities 
f(NA,Ns) from our models. — (a) an irregular 
mass distribution of "cloudlets" (no dependence on 
d), — (b) and (c) the d = 0.01 and 0.1 cases for 
a smooth disk/slab-like structure with column dis- 
tribution Nn(r) oc r -4 , — (d) the d = 0.25 case 
for the smooth spherical structure, — (e) and (f) the 
d = 0.5 and 1.0 cases for the smooth disk/slab-like 
structure. For the smooth mass distribution models, 
note the increased probability density "ridges" astride 
the Na = Nb correlation line, and that as d increases, 
the (Na,Nb) pairs become less correlated with the 
maximum column density of a correlated pair given 
by N H (R)(d/2)-\ 
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impacting LOS is likely to have small Nh (hit at large 
r on the disk). At d > 1 the /(# ac ) oc AT" 1 - 5 , and 
are indistinguishable from the N ac distribution arising 
from irregular structures. 

In Fig. 3, we show the distribution of velocity dif- 
ferences f(AV = |Va — Vb\) between LOS for the 
isotropic sphere (d independent), the rotating disk, 
and the planar inflow/outflow slab. If a significant 
systematic component to the velocities exists in the 
disk/slab, an observable change in /(AV) will oc- 
cur as a function of d, with /(AV) highly peaked 
around AV = for small d. The distribution shapes 
of the rotation and planar inflow/outflow cases are 
virtually indistinguishable for d <J 1, except that 
the maximum velocity difference for the planar in- 
flow/outflow case has an inclination dependent cutoff 

at AV max = 2V* {l — (d/2) 2 } 1 2 . As d increases over 
the range 1 <J d < 2, the decreasing AV max for the 
planar inflow/outflow case becomes clearly distinct 
from the rotation case with AV max = 214. 

Diagnostics from DLOS Observations 

As an example of an application of the tests in 
the previous section, we consider the current obser- 
vational data at z ~ 2 (noting that these data are 
limited by small sample sizes, confusion due to blend- 
ing, and uncertain Nh on the flat part of the curve of 
growth). Our maximum likelihood analysis of the D94 
data (following McGill (1990)) yields a most proba- 
ble size of R = 74 h~ x kpc for spherical absorbers, 
and R = 106 kpc for disks. Adopting this value 
for the size of our single population of absorbers, we 
have < d < 0.025 over the observed redshift range 
(2.1 £ z £ 2.7) along UM673 A,B, and at z ~ 2, we 
have d ~ 0.1 for Q2345+007 A,B 6 and d ~ 0.5 for 
Q1343+266 A,B. 

From coincident lines, the /(NajNb) distinguish 
between structures with smoothly varying monoton- 
ically decreasing Njj(r) oc r -4 and those with irreg- 
ular Nh distributions. Figure 1 (b-e) shows the pre- 
dicted progression from strongly correlated /(Na , Nb) 
for UM673 A,B to relatively scattered for Q1343+266 
A,B, which is generally consistent with the observa- 
tions (F84, S92, B94, D94). The level of f(N A ,N B ) 
correlation in the UM673 A,B data (Fig. 5a of S92) 



6 Fischer et al (1994) and Steidel&; Sargent (1991) find evidence 
that Q2345+007 A,B is a lensed pair. If it is not lensed, then 
we have d ~ 0.5, similar to Q1343+266 A,B. 
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Fig. 2. — The distribution of anticoincident column 
densities f(N ac ). All distributions have been normal- 
ized at Nh(R)- Each curve is labeled by its corre- 
sponding d parameter. — (a) the smooth spherical 
absorbing structures. — (b) the smooth disk/slab- 
like absorbing structures. 
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Fig. 3. — The distribution of LOS velocity differences 
f(AV = \Va — Vb\)- All distributions have been nor- 
malized at AV = 0. Each curve is labeled by its 
corresponding d parameter. The solid curves are for 
smooth disk/slab-like structures with constant rota- 
tional velocities. The dashed line is for disk/slabs 
with constant planar inflow/outflow velocities. The 
distribution arising from planar inflow/outflow ex- 
hibits a sharp cut off at AV/V* = 2 {l - (rf/2) 2 } 1/2 , 
and becomes distinguishable from the rotating disk 
for d > 1.0. The dotted line is for smooth spherical 
structures with isotropic "cloudlet" velocities, which 
give d independent distributions. 
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is not consistent with Nh distributed irregularly on 
kpc scales. The Q2345+007 A,B data (Fig. 6 of S92) 
and the Q1343+266 A,B data (Fig. 3 of D94) are 
scattered about the correlation line, consistent with 
smoothly varying Nh distributions on scales of tens of 
kpc within Lya absorbing structures, and not with an 
irregular distribution. However, a higher level of cor- 
relation for large d cases could arise if we relax the as- 
sumption of a single population such that each struc- 
ture has an Njj gradient smaller than r -4 . Scatter 
would be added to model predictions if we included a 
range of d values resulting from variation in structure 
size and in DLOS separation with redshift. 

From anticoincident lines, the f(N ac ) distinguish 
between irregular, smooth spherical, and smooth disk/ 
slab-like structures for d < 1. As seen in Fig. 2, 
smooth spheres yield a vanishingly small fraction of 
anticoincidences at large column densities, while for 
d <J 0.2 smooth disks yield a significantly smaller frac- 
tion of anticoincidences than do irregular Nh distri- 
butions. For example, if d = 0.1 we find ~ 12% of the 
anticoincident lines that arise from smooth disks have 
Nh > 2Nh(R), whereas ~ 71% are above this Nh in 
the case of irregular distributions. Significant differ- 
ences between smooth spheres and disks persist for 
larger d values. For d = 0.5, no anticoincidences with 
N ac greater than ~ 15Nh(R) should be observed if 
the geometry is spherical as compared to 26% for disk 
models. Since, for d = 0.5, ps 44% of the absorption 
lines are anticoincidences, ps 11% of the total num- 
ber of absorption lines should have N ac > 15Nh(R) 
in the case of a disk geometry. For small samples, 
DLOS separations in the range 0.25 < d < 0.5 are 
likely to provide the greatest leverage. 

In current studies, there are two anticoincident 
columns from which we infer that a single population 
of smooth spherical Nn(r) Lya absorbers is not con- 
sistent with the data. First, it is important to account 
for the fact that in an equivalent width limited sample 
the minimum detectable Nh varies with wavelength 
(redshift). The minimum requirement for an anticoin- 
cidence is that a line of equal strength can be detected 
in the other LOS within a defined velocity/redshift 
window at the level chosen for the analysis. But even 
if the Nh detection limit in a LOS is an order of mag- 
nitude below a detection in the other LOS, a biasing 
toward anticoincidences with large column densities 
may be occurring. For the d = 0.1 Q2345+007 A,B, 
F84 find an anticoincident line in spectrum A (line 
Al) with N H = 3 x 10 15 cm" 2 (6 = 30 km s" 1 ). The 



B spectrum limiting column in this velocity/redshift 
window is Nh = 6 x 10 14 cm -2 , a factor of 5 lower. 
With D ~ 10 kpc, a smooth sphere with Nn(r) oc r -4 
should have a coincident line in spectrum B in the 
range 10 15 < Nh < 10 16 cm -2 . Here is a clear case 
in which an anticoincident line of N ac / Nh(R) ~ 5 is 
seen when the f(N ac ) cut-off is N ac / Nh(R) ~ 1.5 
(see Fig. 2). For this anticoincident line to be con- 
sistent with a single population of spheres having 
Nn(r) oc r~ a , the column distribution would have 
to fall unreasonably steeply, with a > 15. A simi- 
lar argument holds for anticoincidence Bl from B94. 
Caution is in order, however, because an analysis of 
anticoincidences that relies on such small numbers is 
suspect to contamination of the Lya sample by metal 
lines (Steidel & Sargent 1991). 

From coincident lines, the f(AV) diagnose the de- 
gree of systematic motion in Lya absorbers through 
analysis of correlations of the RMS velocity differ- 
ences (AV TTns ) with d. Both the UM673 A,B (d ~ 
0.01) and Q2345+007 A,B (d ~ 0.1) data yield AV rms 
consistent with km s" 1 (S92, F84). The ob served 
AV TTns of 17 km s" 1 (Fig. 10 in S92) for UM673 A,B 
implies that V* <^ 10 km s _1 , unrealistically quiescent 
for an isotropic sphere. However, a 10-20 km s _1 
spread is plausible for disk/slab-like structures in 
which the DLOS penetrate through a flattened ge- 
ometry at some inclination. Systematic motion over 
large scales in the plane should result in a correla- 
tion of AV Tms with d, but due to the 25 km s _1 er- 
ror in the Q2345+007 A,B data, it is not possible to 
determine if the RMS velocity spread observed for 
Q2345+007 A,B is larger than that for UM673 A,B. 
However, the larger Al/ rms = 65 km s _1 observed for 
Q1343+266 A,B (D94) does provide evidence for sys- 
tematic motion. It should be possible to distinguish if 
the motion is predominantly rotational or planar in- 
flow/outflow from DLOS with d ;> 1.5, since the latter 
should exhibit a maximum cut off Al/ max < 214- 

Single LOS data can be compared to DLOS data 
to obtain further diagnostic information. Cowie et al. 
(1995) have resolved Civ AA1548, 1550 lines associ- 
ated with Lya forest absorbers using HIRES (Vogt 
et al. 1994) on the Keck 10-m telescope, finding that 
the absorbers are comprised of a few cloudlets along 
the LOS. The cloudlet velocity dispersions a\ os ps 
18 km s _1 along a single LOS are notably less than 
A 14ms = 65 km s _1 , measured for the D ~ 40 kpc 
pair Q1343+266 A,B (D94). The inequality AV rms > 
cios is not consistent with a spherical absorbing struc- 
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ture in which multiple cloudlets are moving isotropi- 
cally (AV TTns of the cloudlets between DLOS average 
out to be less than a\ os for the cloudlets along a single 
LOS). Many other models, such as spherically sym- 
metric infall or an asymmetric velocity distribution 
for a spherical distribution of cloudlets, are not fully 
consistent with the data. The observed AT/ rms > a\ os 
does, however, arise naturally in a model in which the 
motion is confined within the plane of a disk/slab-like 
structure. 

The detailed distribution in position and velocity 
of the cloudlets responsible for Lya absorption can- 
not be uniquely determined on the basis of current 
observations. However, some models are ruled out 
by our tests, while others remain consistent with the 
data. In particular, DLOS passing through cloudlets 
(blended in Lya) embedded in randomly inclined ab- 
sorbing structures with flattened geometries and a 
smoothly declining Nn(r) can give rise to the ob- 
served /(NajNb) and f(N ac ) distributions, and to a 
LOS cloudlet velocity dispersion a\ os = 18 km s _1 . If 
the kinematics is systematic as we have modeled, then 
the observed correlation (tentative) of AT/ rms with d 
will hold true and the f(AV) will generally behave as 
shown in Fig. 3. 

Conclusions 

The new generation of 8-m class telescopes and 
efficient spectrographs will result in a more than ten- 
fold increase in the number of coincident and an- 
ticoincident lines, column densities will be directly 
measured, and blends will be separated. Expanded 
analyses similar to those discussed in §3 can then 
be pursued, comparing model predictions to the col- 
umn distribution of pairs, /(NajNb) (Fig. 1), the 
distribution of anticoincident lines, f(N ac ) (Fig. 2), 
and the distribution of velocity differences between 
pairs, f(AV) (Fig. 3). The discovery of more quasar 
pairs, and even groups, with separations ranging from 
rf=0.01toc?=1.5, will allow the full range of diag- 
nostics to be applied. 

The simple models that we have considered here, 
though idealistic, should bracket the range of rea- 
sonable conditions comprising Lya forest structures. 
These models provide clear predictions for correla- 
tions between the d parameter and the distributions 
of coincident and anticoincident column densities and 
velocity differences. Based upon a loose interpreta- 
tion of current data, a flattened structure with a com- 



ponent of systematic velocity (either rotation or pla- 
nar inflow/outflow) is favored over a spherical struc- 
ture with isotropic velocities. The large differences 
in the predicted distributions for the extreme models 
provides great optimism that the elusive answer to 
the questions of taxonomy, cosmogenesis, and general 
nature of Lya forest clouds may soon be within our 
grasp. 
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